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Abstract 



The renormalized coupling g R defined through the connected 4-point function at 
zero external momentum in the nondinear 0(3) sigma-model in two dimensions, 
is computed in the continuum form factor bootstrap approach with estimated error 
~ 0.3%. New high precision data are presented for g R in the lattice regularized theory 
with standard action for nearly thermodynamic lattices L/£ ~ 7 and correlation 
lengths £ up to ~ 122 and with the fixed point action for correlation lengths up to 
~ 12. The agreement between the form factor and lattice results is within ~ 1%. We 
also recompute the phase shifts at low energy by measuring the two-particle energies 
at finite volume, a task which was previously performed by Liischer and Wolff using 
the standard action, but this time using the fixed point action. Excellent agreement 
with the Zamolodchikov S-matrix is found. 



1 On leave from the Institute of Theoretical Physics, Eotvos University, Budapest, Hungary 



1 Introduction 



The presently only known practical way to define a relativistic quantum field theory 
non-perturbatively in 4 dimensions is by using the lattice regularization. For example 
it is hoped that one will, once sufficiently powerful computers are available, be able 
to answer the question whether QCD is the correct theory of the strong interactions 
by studying the continuum limit of the lattice theory. 

It is however notoriously difficult to control the continuum limit of a lattice 
theory. Analyses of lattice Feynman graphs, as initially performed by Symanzik ]IJ], 
show that order by order in renormalized perturbation theory, physical quantities 
approach their continuum limit as integer powers in the lattice spacing (up to loga- 
rithmic corrections). However, since it is not known that such behavior controls the 
approach to the continuum limit in the full non-perturbatively defined theory, the 
invocation of such a power-law approach to extrapolate data produced in numerical 
simulation experiments has merely the status of a (plausible) working hypothesis. 
The integer nature of the powers adopted in studies of theories such as QCD, is 
considered to be connected to the widely expected property of asymptotic freedom. 
Here again, the very question whether the continuum limit of the lattice theory 
really describes an asymptotically free theory is highly non-trivial.0 

This work is part of an on-going effort of the present authors to test whether 
the "conventional wisdom" is correct in a simpler model, the non- linear 0(3) sigma 
model in 2 dimensions. This model is, like QCD, perturbatively asymptotically free 
and also has instanton and superinstanton solutions. It has however classically 
the additional beautiful property of being integrable; in particular there exist an 
infinite set of non-local conserved charges. Assuming that the quantum theory has 
a mass gap and the spectrum contains a vector multiplet of stable particles, the 
existence of such conserved charges in the quantum theory forbid particle production 
in scattering and, as shown by Liischer || fixes the 2-particle S-matrix to that 
previously postulated by Zamolodchikov and Zamolodchikov (ZZ) [|j (up to CDD 
ambiguities) . 

All these properties were obtained starting from a formal Lagrangian, where one 
first computes off-shell Appoint functions and goes on-shell via the LSZ formalism to 
obtain the S-matrix elements. The so-called form factor bootstrap (FFB) approach 
HIJIQ P rocee ds in the other direction. One attempts to obtain off-shell information 
starting from the knowledge (postulate) of the stable particle spectrum and their 
S-matrix. In the first step one constructs the form factors of (composite) operators, 
satisfying all physical constraints (analyticity, generalized Watson theorem etc.) and 
then Green functions are obtained by saturating with a complete set of states. This 

2 Indeed the authors of this paper are divided into two subsets having different opinions on the 
probable answer! 
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is a program which is only feasible in a theory where there is no particle production 
(i.e. only in 2 dimensions). Even this program, for which there has recently been a 
lot of progress [|qi9yi0| , involves mammoth effort. 

Unfortunately since lattice regularization breaks these conservation laws and 
the FFB relies on some nontrivial assumptions, the expectation that the continuum 
limit of the lattice 0(3) model coincides with the FFB is not guaranteed. A first 
investigation of this issue was by Liischer and Wolff [jll]] who computed the phase 
shifts on the lattice by measuring two particle state energies in finite volume. Their 
results (taking account of lattice artifacts) were completely consistent with the ZZ 
S-matrix. In the course of a similar investigation of the nature of the continuum limit 
of the 0(2) model |l^], in order to test our programs and a slightly modified form of 
the analysis, we repeated the measurement of the phase shifts in the 0(3) model. We 
also made simulations with a fixed point action [13|, and found indications that the 
lattice artifacts are much smaller than in the case of the standard action. An account 
of our investigations is given in sect. 4. Our results are again in good agreement 
with the ZZ S-matrix Q. 

An off-shell quantity of physical interest is the current-current ( J-J) correlation 
function. It has been shown that J-J computed in the FFB approach ]^||| agrees 
well with conventional renormalized perturbation theory at high energies, at least 
up to p/M ~ 10000, when the analytical value of the ratio of the mass M to 
the A-parameter [14] is used. Connected to this are two important inter-related 
properties; firstly the ZZ S-matrix shows an "on-shell form of AF", in the sense 
that the phase shifts fall logarithmically to zero with the energy at high energy (see 
sect. 4). Secondly the thermodynamic Bethe ansatz, which is used to compute M/A, 
reproduces the 2-loop (3- function coefficient. The agreement between the FFB and 
lattice computations of J-J is also within ~ 1% for the entire range of momenta up 
to ~ 40M. 

Despite this wealth of circumstantial evidence for the validity of the conventional 
picture, there is still room for doubt. In particular in a recent paper, two of the 
present authors [15], assuming a certain form of the lattice artifacts, have found 
statistically significant deviations between the continuum limit of the lattice J-J 
correlation function and the FFB at low energies. 

Unfortunately the J-J correlation function at low energies is a quantity which 
behaves qualitatively similarly to that in a free theory. It is plausible that a dif- 
ference between two theories would manifest itself more clearly in a quantity which 
vanishes in the free theory e.g. the zero momentum coupling g R defined through 
the connected 4-point function. There is an enormous literature on the computa- 
tion of this quantity in the 2 (and higher) dimensional non-linear sigma models, see 
e.g. ref. |l(| and references therein. The main new contribution of this paper, the 
computation of g R by the FFB, is outlined in sect. 2. This is the first time that this 
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method has been used to compute a 4-point function, and it is rather surprising that 
one is apparently able to get such a good approximation for g R . 

In sect. 3 we present results on g K using two different lattice regularizations, 
the standard action (including new high precision data on thermodynamic lattices 
at large correlation lengths) and the fixed point (FP) action. The nature of the 
approach to the continuum limit is not so clear, but whatever (reasonable) extrap- 
olation is made, it agrees with the truncated FFB result better than ~ 1%. 

2 Computation of g R in the continuum theory 

There have been various approximation schemes to compute low energy (non-pertur- 
bative quantities) in the continuum formulation of the 0(n) models in two dimen- 
sions: the (/-expansion jl7]], the e-expansion [jl^], and the \/n expansion pOj . In 
this section we will present a new approximation using the form factor bootstrap. 

2.1 Definitions 

Consider a general continuum quantum field theory in two dimensions, in infinite 
volume, with global O(n) symmetry. Let a a (x), a = 1 ... n be a vector multiplet of 
(renormalized) Euclidean scalar fields with 2-point function 

S a ^ 2 (x 1 ,x 2 ) = (a a ^x 1 )a a2 (x 2 )). (2.1) 
The inverse of its Fourier transform 

G(k)5 aia2 = j d 2 xe ikx S aia2 (x,0), (2.2) 
is assumed to have an expansion for small momenta 

G{k)- 1 = Z-^Ml + k 2 + 0(k 4 )^j . (2.3) 
We denote the 4-point function by 

S a ^ a ^( Xl ,x 2 ,x 3 ,x 4 ) = <^iK 2 (x 2 K 3 (* 3 K 4 (*4)}, (2.4) 
and the connected 4-point function 

S2 ia2a3a4 (x 1 ,x 2 ,x 3 ,x 4 :) = 5 aia2a3a4 (xi,x 2 ,X3,x 4 ) - 5 aia2 (xi,x 2 )5 a3a4 (a;3,X4) 

-S a ^( Xl ,x 3 )S a2a4 (x 2 ,X4) - S aia4 (x 1 ,x 4 )S a2a3 (x 2 ,x 3 ). (2.5) 

Introducing the Fourier transform by 

, 4 

S a ^ a ^(h,k 2 ,k 3 ,h) = U [d 2 x je ik ^] S^^ixux^x^xt) , (2.6) 
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and similarly S c for the connected part, the conventional zero momentum (dimen- 
sionless) coupling (in two dimensions) is defined by 



furl i 



G(0) 2 n 2 

a.b 



where 



££i«a«3«4( fclj fe2j fc3> fc4 ) = (27r) 2 5( 2 )(A;i + fc 2 + A; 3 + W 11 ™ 4 ^!, fc 2 , fc 3 , fe 4 ) . 

(2.8) 

We will assume that the 2-point function has a spectral representation 

° {k) = z r ' (2 - 9) 

where the normalization constant Z takes into account that, assuming that the 
spectrum of the theory contains a vector multiplet of stable particles of mass M, we 
normalize the spectral density p so that the 1-particle contribution is 

pW(p) = S( f i-M). (2.10) 

Then the coefficients appearing in the small momenta expansion above can be ex- 
pressed as 

Z R = Z-f, (2.11) 
o< 



2 



M 2 6 2 



(2.12) 



where 72 and <5 2 are the moments 

72 =M 2 |d^, (2.13) 



h -V J / dM#. (2.14) 



Further the coupling g R can be written as 

+ 2 

n 72 (5; 

where 74 is defined through 



= T~ (2-15) 



V-2 

G^ * 4 * * (0,0, 0,0) = —J± ^s a ^S a3a4 + S aia3 S a2a4 + S aia4 S a2a3 ) . (2.16) 
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Strictly speaking, all observable physics (in a massive theory) is on-shell; dif- 
ferent interpolating fields give the same results. Off-shell amplitudes of particular 
composite (and 'elementary') operators, are only of physical interest if they are 
sources of (idealized) infinitesimal weakly interacting probes. The fields are charac- 
terized by their various quantum numbers and their dimensions e.g. coded in the 
behavior of the 2-point function G(k). The assumption eq. (|2.9|) corresponds to a 
limitation on the behavior of the associated spectral function p(/j,) as fi — * oo. It 
is an implicit connection to the association of a a {x) with a particular local field 
in the Lagrangian quantum field theory. In the non- linear 0(n) sigma model the 
'elementary' field a a (x) is characterized by its vanishing engineering dimension (in 
addition to being an isovector and spacetime scalar). We can therefore assume that 
the short distance singularities of its two-point function are sufficiently weak (only 
logarithmic) so that the spectral representation ( [2.9| ) holds without subtractions. 
Indeed, in the FFB approach, this uniquely defines (up to normalization) an oper- 
ator <7 a (x), which has form factors that are not growing too fast at infinity so that 
the corresponding spectral density p(/i) vanishes fast enough for jj, — > oo. 

2.2 1/n- and e-Expansions 

The leading order computations for the spectral integrals in the 1/n-expansion have 



been performed in ref [19 



72 = 1 + 0.00671941- +0 -7 , (2.17) 



1 


+ o( 






v n 2 ) 


n 






1 


+ 01 









S 2 = 1 + 0.00026836- + O , (2.18) 

n In 2 / 



and also for the coupling [ 20 1 



g R = — ( 1 - 0.602033^ + 0(4r) ) , (2.19) 



n \ n n 2 



which gives the approximation g R = 6.70 for the case n = 3. 

In the ^-expansion one obtains g R = 6.66(6) []17| and the e-expansion g R = 
6.55(8) |j"6]| . Considering the rather short series in each case it is amazing how well 
the estimates by the various methods agree. 



2.3 Form factor bootstrap computation for n = 3 

The form factor bootstrap aims at reconstructing iV-point functions of local opera- 
tors of integrable field theories from the knowledge of the spectrum of stable particles 
and their S-matrix. A description can be found in Smirnov's book ||, the review of 
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Karowski ||, a recent paper [ 10 1 and in various articles of two of the present authors 

To our knowledge this is the first time that the method has been applied to the 
computation of 4-point functions. The computation is rather involved and here we 
will only give a very brief outline and present our results. The calculation in this and 



in other integrable models will be described in detail in a forthcoming paper [24]. 

We assume (as did the Zamolodchikov brothers in their construction of the 
S-matrix) that there are no bound states in the O(n) models. Then the spectral 
density p has an expansion over contributions from the intermediate states with an 
odd number of particles (due to internal parity symmetry) 

oo 

P {n) = Y;p {2k+1) ^)> ( 2 - 2 °) 

k=0 

and correspondingly the spectral integrals 

oo 

(2fc+l) 



72 = 1 + I>r +1) > (2.21) 

k=l 

oo 

5 2 = 1 + J2^2 k+1) ■ (2-22) 



fc=i 



The form factors are given by Smirnov for the 0(3) model and have been 
recomputed in ||. For example the matrix element of the Minkowski operator <r a (0) 
associated to the Euclidean field cr a , connecting the vacuum to the 3-particle in-state 
with rapidities 0\ > 2 > 3 , is given by 



with 



(O\a a (O)\a 1 ,0 1 -a 2 ,0 2 ;a 3 ,e 3 r = Vz^ ia2a3 (0 u 2 ,0 3 ) (2.23) 

^ ia2a3 (0i,02,0 3 )=K 3 n i<] m-O j ) 
h ~ 2 )5 a ai 5 a2a3 + (0 1 - 3 - 27ri)S a a2 6 a3ai + (0 2 - 5 ai J , (2.24) 



where 



The expression of the 5-particle matrix element is also known explicitly but it is much 
too long to be written here. Using these results the 3- and 5-particle contributions 



6 



r 


r 

72 


6 (r) 
°2 


3 


1.67995(1) x 10" 3 


3.46494(1) x 10" 5 


5 


6.622(1) x 10" 6 


7.114(1) x 10~ 9 



Table 1: r-particles contribution to 72,^2 

to 72 and 5 2 are given in table |. It seems that the series converge extremely rapidly 
and we would estimate 

72 = 1.001687(1) , (2.26) 
5 2 = 1.000034657(1) , (2.27) 

where the estimated errors come from inspecting the pattern of relative n-particle 
contributions suggested by the 1,3,5-particle states. 

The 4-point function has an expansion in terms of contributions of intermediate 
states with l,m,n particles respectively 

S aiam (k u k 2 ,k 3 ,k 4 ) = (2TT) 2 6^(k 1 + k 2 + k 3 + k 4 ) 

■ M -6 y apiap2ap3ap4 (kpi,kp 2 ,kp 3 ,kp A ) (2.28) 

perms P 

where h = (k a , k i2 ), 

V^ {kl MMM) = {2,)Hf £ m + ^)^ + h 1 + k21 ) 

f— 1 Ei - ik u E m - tk u - ik 22 

• T\~^ 0|^H0)I1)(^ , (2.29) 

En + ik 42 

where n run over all states with odd numbers /, n of particles, and m over all states 
with an even number m > 0. The somewhat symbolic ^ in ( [2.29D really means in 
addition to the summation over all internal quantum numbers of the particles and 
integration over all particle rapidities a sum over the integers I, m and n. The limit 
of zero momenta is very delicate because each term in the sum is a distribution in 
the momenta where the singularities occur when certain linear combinations of the 
momenta are zero. In particular the contributions from terms in the above sum with 
m = 0, not only cancel the disconnected pieces S — S c , but also produce extra terms 
proportional to 8(kn + k 2 i). The singularities must be canceled by other terms in 
the sum with m > e.g. the singularity from the contribution 1-0-1 is canceled by 
that of the contribution 1-2-1. We can avoid this problem if, as is assumed in the 
following, we restrict ourselves to momenta where kn + kji 7^ for any i ^ j. An 
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1, m, n 


T4;2mn 


1,2,1 


-4.16835(1) 


1,2,3 


0.05175(1) 


3,2,1 


0.05175(1) 


1,4,1 


-0.004065(1) 



Table 2: Z-m-n-particles contribution to 74 

additional technical complication comes from the fact that each term in the sum is 
rather involved because the matrix elements have parts with differing connectivity 
properties, e.g. the matrix element occurring in the 1-2-1 contribution (9\ > 62)'- 

Z-5(b, 0|^(O) |a a , #1; as, 02) [n = n axa ^ + i7r- ie, 9 1 ,9 2 ) 

+4n5^6 bai 6^ - 0!) + Att5(<P - 9 2 )S aiariab (0 ~ O2) , (2.30) 

where the S-matrix elements S ai a 2 ;ab are given in sect. 4. 

The practicability of the computation of the zero-momentum coupling using the 
form factor bootstrap approach obviously depends crucially on the question whether 
the sum over intermediate states for the 4-point function, 

74 = ^2 ^;lmn (2.31) 
l,m,n 

converges rapidly. We startedP] to investigate this question in the Ising model and we 
found that the I + m + n = 6 contributions are much smaller than the leading 1-2-3 
term pjfl . Fortunately, for the 0(3) model under investigation here the situation is 
rather similar. The contributions of the l-m-n intermediate states with l + m + n < 6 
to 74 are given in table [2|. 

The leading 1-2-1 contribution is a factor ~ 42 greater in magnitude than the 
sum of l-m-n contributions with l-\-m + n = 6. It is extremely difficult to bound the 
rest of the contributions, especially since the signs are not known in general. Even 
the computation of the states with / + m + n = 8 would be quite an undertaking. 
But assuming that the pattern in table ^ continues (as for the case of the 2-point 
function) and that the sum of the remaining contributions I + m + n > 8 is < 10% 
of the sum of the / + m + n = 6 contributions, we obtain the result 

74 = -4.069(10) , (2.32) 

and hence our final result 

g R = 6.770(17) . (2.33) 
3 We hope to be able to work out the I + m + n = 8 contributions in the Ising model soon. 
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3 Lattice computations of g R 

In the framework of the lattice regularization there are two methods to compute g R 
in the O(n) models. The first is using high temperature (strong coupling) expansions 
and the second by numerical simulations. 

For numerical simulations we consider a square lattice with both the standard 
action 

S = -(3^2a(x)-a(x + jl), (3.1) 

X,fi 



where cr(x) ■ cr(x) = ^2 a o' a (x)a a (x) = 1 and the fixed point (FP) action of ref. 13] 



3.1 High temperature expansion 

Concerning the high temperature (HT) expansion for the standard action, long series 
flg[| have been obtained for the 2- and 4-point functions and for the second moment 
H2 defined by 

M2 ^ ia2 = x 2 (a ai {x)a a2 (0)) (3.2) 



to obtain M R through M R = 4G(0)//i2, where G(k) is defined analogously to eq. (2.2) 
but with the integral replaced by a sum. The analyses of the HT expansion for the 
spectral moments give 72 = 1.0013(2) |J| and 5 2 = 1.000029(5) |l||. The agreement 
with the FFB values eqs.( |2.26| , |2.27| ) is acceptable; note that these are smaller than 
that anticipated from the leading order of the 1/n approximation, eqs. ( |2.17 ), ( 2.18 ). 
The coupling is defined as the continuum limit 

g R = lim g R (£) , (3.3) 

where £ is the correlation length in lattice units. The analysis is hampered by 
the lack of rigorous knowledge of the position of the critical point and the exact 
approach to it. In particular, the conventional wisdom that the critical point is at 
(3 = 00 is usually built into the analyses. The various Pade approximations show 
the coupling falling rapidly as (3 increases in the region of small /3, then a region of 
rather flat behavior after which the various approximations show diverse behavior 
(see e.g. fig. 12 in ref. (2^]) making error estimates rather difficult 



In ref. [23] Campostrini et al quote for the case n = 3 the result g R = 6.6(1), 

and in a more recent publication Pelissetto and Vicari cite 6.56(4) |l6[ . Butera and 

Comi on the other hand are rather cautious, and did not quote a value for the case 

n = 3 in ref. [T^j; if pressed they would at present cite g n = 6.6(2) p5[ . 

4 A result is considered more reliable if in the region where the coupling flattens the approximants 
have no complex singularities near the real axis. Unfortunately for the special case n — 3 there do 
tend to be nearby singularities |js||. 
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3.2 Numerical simulations 

Monte Carlo computations of g K have a long history, see e.g. refs. |L7],^(|. In order to 
attempt to match the apparent precision attained in the FFB approach outlined in 
section 4, we decided to perform even more precise measurements than were carried 
out previously. 

We work here on a square lattice of size L in each direction and periodic bound- 
ary conditions. The coupling in finite volume is defined through Binder's cumulant 

Sa(0= lim 9k(£,L), (3.4) 



2 

1 + - 

n 



2\2\ 



(J?) 



2\2 



(3.5) 



where S a = ^2 X a a {x). In this definition we have taken (as in ref. |T 



i 



/ G(0) 

2sin(^/L)\/ G(k ) 



(3.6) 



where &q = (2w/L,0). In this paper we will use £ to denote the second moment 
correlation length, to which ( |3.6[) converges, for large L. Although conceptually 
different, £ is very close to the exponential correlation length (£ ex p): also using (2.3) 
and ( [242] ) in the definition the FFB results pT26[ ) and ( pT27j ) yield 



e 



cxp 



'21 = 1.000826(1) 

02 



(3.7) 



The standard action Monte Carlo measurements were performed using a method 
similar to the cluster estimator of [11|. We measured g R at correlation lengths 
ranging from 11 to 122 at L/£ ~ 7. These measurements were used to investigate 
the approach to the continuum. To study the finite volume effects, we repeated the 
runs at £ ~ 11 on three other lattices with L/£ ~ 5.5, L/£ ~ 9 and L/£ ~ 13. 
The results of all these runs (together with the preliminary results corresponding to 
£ ~ 167, 230 and 309) are recorded in table § 

In this table we also indicate the number of measurements. Each run consisted 
of 20k sweeps of the lattice with measurements after each sweep. The error was 
computed from this sample by using the jack-knife method. 

We have measured g n with the FP action at three different values of 0: 0.70, 
0.85 and 1.00, corresponding to correlation length £ ~ 3.2, 6.0 and 12.2, at the 
values of z = L/£ in the range 5.4 . . . 8.2. 
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To get a feeling of the finite volume effects we took the expression for g R in 



the leading order 1/n expansion from [2C] and simply replaced the integral over 
momenta by a discrete sum. In this way we obtained 



g R [L] = <7 R [oo][l - aiv/Ie^(l + 0(1/ z)) + 



(3.8) 



with a\ = \/8n = 5.013 for large z = L/£. 

Figure || shows these results for g R plotted against the combination sfz exp(— z), 
motivated by the 1/n result fl3.8| ) . 

We have determined the MC prediction of g R both for the standard action and 
the FP action. Making a linear fit in y/zexp(— z) to the four data points at (3 = 1.50 
for the standard action (see fig. [l]) one obtains 



g R (z = oo,£ = 11.0) =6.60(1) 



(3.9) 



with a\ = 5.0(4). We therefore used ( |3.8p to "renormalize" the data with z ~ 7 in 
table H to the common physical size z = 7.25. 

Next we extrapolated the results of our measurements to the continuum limit. 
In Fig. H we show both the measured and the corrected values of g R versus l/£ for the 
data with L/£ ~ 7. It clearly indicates that the continuum limit of g R is approached 
from below as is the case indicated by some Pade analyses and in the leading order of 
the 1 /n approximation of the lattice theory |20] . The only theoretic framework for 
estimating lattice artifacts comes from considering Symanzik's j[| effective action. 
The absence of even parity 0(3) invariants with odd engineering dimensions suggests 
an approach to the continuum limit with leading behavior (log£) r /£ 2 . As stressed 
in the introduction, there is no rigorous non-perturbative proof of this behavior. 

Figure |2| shows two fits. The first is a quadratic fit of the form suggested by 
the Symanzik analysis, where we have taken r = 1 0: 



5r (£, z = 7.25) =g R (oo,z = 7.25) 



, h logg 62 

+ £2 + £2 



(3.10) 



with g R (oo, z = 7.25) = 6.702(16), b x 
fit in the fig. § of the form: 



-4.4(2.2), b 2 = 8(5). We made also a second 



g R (£,z = 7.25) =g R (oo,z = 7.25) 



1 + 



(3.11) 



with g R (oo,z = 7.25) = 6.710(13), d 2 = -0.27(4). Although there is no theoretical 
basis for such a fit, it describes the present data as well as the first. 



5 this is the analytic form found in the the leading order 1/n expansion with b 1 
-|ln2 + i @ 



b> 
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The morale is that the precision of our measurements (which are the presently 
best available) does not really discriminate between the two fits (or between any 
intermediate fits with leading behavior (ln£) r /£ 2 with say large r). This unfortu- 
nately limits the accuracy on the MC determination of the continuum value of g R ; 
a conservative estimate would be g R (oo, z = 7.25) = 6.71(2). 

6.9 



6.8 



6.7 



6.6 



6.5 



6.4 



6.3 



6.2 





P FP =0.70 




O P FP =0.85 




□Ppp=1.00 




■ P ST =1.50 




P ST =1.60 




AP ST =1.70 




P ST =1.80 




• P ST =1.90 


*FFB 






1 


1 \ , 



0.005 



0.01 



0.015 



1/2 , . 

z exp(-z) 



Figure 1: The coupling gn(z,[3) for the FP (open symbols) and standard actions (filled symbols) 
plotted vs. y/zexp( — z). The correlation length for the FP action are in the range 3.2 . . . 12.2 while 
for the standard action in 11 . . . 122. The theoretical value from the FFB is also shown. The linear 



fits are motivated by the the 1/n prediction (3.8) and are based on the (3st = 1-50 and /3fp = 0.85 
data, respectively. 



Finally, we used ( |3.8| ) again to extrapolate the continuum result for the finite 
physical size z = 7.25 to the thermodynamic limit z = oo. Thus our final result 
based on the standard action lattice measurements is: 



6.77(2) . 



(3.12) 



This error contains both the ambiguity in using the fits ( |3.10 ) or ( |3.11| ) and the error 
in oi< 

The data for the FP action seem to lie on a universal curve (the slope of which 
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□ original data 

• corrected data 
1 fit 

x extrapolation 
1 fit 

O extrapolation 




Figure 2: The measured values for £tr together with the correc ted o nes at L/£ = 7.25 for the 
standard action. The solid and dotted lines correspond to the fits (3.10) and (3.11) respectively. 



roughly corresponds to the 1 jn prediction) in spite of the extremely short correlation 
length. We interpret this as an indication that the lattice artifacts for the FP action 
are small, in any case smaller than our error bars. The measured values for the 
standard action show a considerable lattice artifact but the data at z ~ 7 seem to 
converge to the FP result with increasing £. Extrapolating the FS effects to z = oo 
we get 



6.77(2) 



(3.13) 



and a\ = 5.0(3). This extrapolation is based on the four /Spp = 0.85 data points but 
including all points of table ^| does not alter the extrapolation significantly. 

Our results ( 3.12) ) and ( |3.13| ) are clearly above the value determined previ- 
ously via finite size scaling |2(|, although statistically compatible with Kim's error 
estimate. They are above the central values quoted on the basis of high tempera- 
ture expansions discussed above but very consistent with each other and the result 
eq.( 2,33| ) from the form factor bootstrap. 

We would like to end this section with a comparison of an analytic prediction 
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1.15 



1.10 



1.05 



1.00- 



I 



1.4 1.6 1.8 2.0 2.2 2.4 

P 



Figure 3: Th e measured values of the x/£ ratio divided by the 4-loop approximation to the 



prediction (3.15) 



of the ratio x/£, 2 with the MC data. In ref. ||] the perturbative short distance 
expansion of the spin 2-point function was refined to 

^-^(^"M^)}- <314) 

The new result was the exact (though non-rigorous) determination of the overall 
nonperturbative constant. Using this we can straightforwardly derive the relation 

X 3tt7| 1 J, , A c 



The first three non-universal perturbative coefficients are known for the standard 
action lattice regularization p?] ]: 

d = 0.1816, c 2 = 0.1330, c 3 = 0.1362. (3.16) 
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p 


L 


runs 




£ 


X 


9n 


1.5 


60 


110x20k 


5.469(10) 


10.97(2) 


173.76(31) 


6.269(20) 


1.5 


80 


344 x20k 


7.253(7) 


11.03(1) 


175.95(11) 


6.553(16) 


1.5 


100 


350 x20k 


9.050(8) 


11.05(1) 


176.51(6) 


6.613(17) 


1.5 


140 


361x20k 


12.68(1) 


11.04(1) 


176.30(5) 


6.560(18) 


1.6 


140 


214x20k 


7.361(15) 


19.02(4) 


447.30(6) 


6.612(17) 


1.7 


250 


367x20k 


7.246(4) 


34.50(2) 


1267.20(57) 


6.665(14) 


1.8 


500 


361x20k 


7.717(4) 


64.79(3) 


3839.07(1.54) 


6.688(15) 


1.9 


910 


101x20k 


7.4395(6) 


122.32(1) 


11884.9(7.0) 


6.738(24) 


1.95 


1230 


12 x 20k 


7.352(20) 


167.30(45) 


20835 (57) 


6.751(78) 


2.0 


1600 


4 x20k 


6.949(25) 


230.26(83) 


36826(142) 


6.981(132) 


2.05 


2100 


5 x20k 


6.804(67) 


308.63(3.06) 


63011(731) 


6.608(216) 



Table 3: f, \ — 3G(0) and gn for the standard action 



Fig. H shows the measured values of the x/£ 2 ratio divided by the 4-loop approxi- 
mation to the prediction (3.15). Note that ( p. 15 ) seems to be satisfied by our data 



to a better accuracy than the analogous one for £ alone [27], but the ratios decrease 
rapidly and we do not know if they eventually overshoot the asymptotic value of 1. 



Note also that the prediction ( |3.15 ) is independent of the M/A^jg ratio 8/e of ref. 
0- 



(3 


L 


runs 






X 


9n 


0.70 


18 


115x360k 


5.682(3) 


3.168(1) 


19.501(8) 


6.493(10) 


0.70 


22 


197 x360k 


6.924(3) 


3.177(1) 


19.646(6) 


6.674(12) 


0.70 


26 


363 x360k 


8.176(3) 


3.180(1) 


19.686(4) 


6.761(12) 


0.85 


32 


31 x 600k 


5.359(5) 


5.971(5) 


55.74(5) 


6.417(14) 


0.85 


34 


189 x360k 


5.678(2) 


5.988(2) 


56.03(2) 


6.485(8) 


0.85 


36 


52 x 360k 


6.006(4) 


5.994(4) 


56.24(3) 


6.581(13) 


0.85 


42 


140 x360k 


6.986(3) 


6.012(3) 


56.51(2) 


6.691(14) 


1.00 


70 


691x100k 


5.734(4) 


12.207(8) 


189.4(1) 


6.487(14) 



Table 4: Data for the FP action. 
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4 Phase shift analysis from 4-spin correlators 

The prediction for the scattering amplitude of two particles at center of mass mo- 
mentum p = Msinh^# in the 0(3) non-linear sigma model by Zamolodchikov and 
Zamolodchikov j|] is given by 

2 

7=0 

where P 1 are the isospin projectors and the phase shifts Si are given simply by 

2i5 (p) _ g + 2 "T , v 

e 2is 1 ( P )_ 6 + 2™ Q-™ U3) 

9 -2m + mt' { ' 



2i5 2 (p) = Z ~ m 

9 + i-K 



(4.4) 



lattice 


P 


Tx L 


m 1 


mL 


A 


1.54 


256 x 128 


13.632(6) 


9.4 


B 


1.40 


128 x 64 


6.883(3) 


9.3 


D 


0.85 


128 x 64 


6.03(1) 


10.0 


E 


0.70 


64 x 32 


3.186(15) 


10.6 



Table 5: Parameters of simulations in the 0(3) model. On lattices A, B the standard action was 
used while D and E denote simulations with the FP action. 



In fact the Zamolodchikovs specified the S-matrix for general n > 3 and these 
expressions have been verified to 0(l/n 2 ) in the 1/n-expansion. For the particular 
case of 0(3), Liischer and Wolff [jllj checked the expressions for low energies using 
MC simulation. The agreement was completely satisfactory; in particular the data 
were consistent with the highly non-trivial non-perturbative property that the S- 
matrix at zero energy is repulsive «Sa'&' ; a&(0) = —^a'b^b'a, which is a crucial condition 
in the FFB construction. We repeated the measurements of [11] for the standard 
action using a modified method of analysis. In addition we also performed measure- 
ments using the FP action. Since the S-matrix is the essential ingredient in the FFB 
construction we report the results here. 

The method of Liischer and Wolff p| is based on the following idea: The 
momentum of one of the particles in a 2-particle state with zero total momentum 
in a periodic box (in 1-dimension) takes discrete values p n given by the periodicity 
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I 


k 


^cxact 


Se 


u VVr 


tn 


M 




o 


1 


1.4595 


1.36(3) 


1.49(2) 


1 


20 


1.36(7) 








1.45(5) 


1.50(4) 


8 


8 






2 


1.2840 


1.25(2) 


1.30(2) 


1 


20 


1.22(3) 








1.38(6) 


1.34(5) 


8 


8 




1 


1 


0.1751 


0.10(1) 


0.17(3) 


1 


20 


0.15(1) 








0.19(3) 


0.19(5) 


8 


8 






2 


0.2692 


0.15(2) 


0.23(2) 


1 


20 


0.23(1) 








0.23(4) 


0.24(5) 


8 


8 




2 


1 


-1.3874 


-1.51(2) 


-1.37(1) 


1 


20 


-1.39(2) 








-1.43(3) 


-1.36(2) 


8 


8 






2 


-1.0944 


-1.11(1) 


-1.06(1) 


1 


20 


-1.05(1) 








-1.04(2) 


-1.06(2) 


8 


8 





Table 6: Phase shifts from lattice A 





k 


<5cxact 


Se 


(5wF 


to 


M 


<>LW 





1 


1.4584 


1.41(1) 


1.48(1) 


1 


20 


1.51(7) 








1.47(1) 


1.49(1) 


3 


10 






2 


1.2818 


1.29(1) 


1.29(1) 


1 


20 


1.2(1) 








1.35(1) 


1.30(1) 


3 


10 




1 


1 


0.1764 


0.09(1) 


0.12(1) 


1 


20 


0.14(1) 








0.12(1) 


0.12(1) 


3 


10 




2 


1 


-1.3859 


-1.42(1) 


-1.35(1) 


1 


20 


-1.38(2) 








-1.36(1) 


-1.35(1) 


3 


10 






2 


-1.0914 


-1.05(1) 


-1.03(1) 


1 


20 


-1.02(1) 








-1.02(1) 


-1.03(1) 


3 


10 






3 


-0.9102 


-0.78(2) 


-0.78(1) 


1 


20 


-0.81(2) 








-0.76(1) 


-0.78(1) 


3 


10 





Table 7: Phase shifts from lattice B 
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I 


k 


^exact 


Se 


5\VF 





1 


1.4726 


1.47(3) 


1.46(1) 




2 


1.3107 


1.29(2) 


1.29(2) 


1 


1 


0.1598 


0.16(1) 


0.14(2) 




2 


0.2545 


0.23(2) 


0.27(1) 


2 


1 


-1.4056 


-1.41(2) 


-1.41(1) 




2 


-1.1321 


-1.15(1) 


-1.16(1) 



Table 8: Phase shifts from lattice D. 



I 


k 


Inexact 


Se 


5\VF 





1 


1.4664 


1.49(1) 


1.48(1) 




2 


1.2977 


1.27(1) 


1.28(1) 


1 


1 


0.1674 


0.18(1) 


0.15(1) 




2 


0.2619 


0.26(1) 


0.28(1) 


2 


1 


-1.3968 


-1.35(1) 


-1.39(1) 




2 


-1.1138 


-1.14(1) 


-1.13(1) 



Table 9: Phase shifts from lattice E. 



condition 

p n L + 25(p n ) = 27m. (4.5) 

Accordingly, the energy of this state is given by 

E n = 2E^( Pn ) = 2^pJTW, (4.6) 

where E^(p) is the energy of a 1-particle state with momentum p. From the 
measurement of the energy spectrum E n for some low lying states one can then 
calculate the momentum p n and using eq. ( [4.5|) the phase shift S(p n ). Varying 
L and taking different values of n one can determine 5{p) at several values of its 
argument. 

To determine the 2-particle energies the correlation matrix has been measured^ 

C xy {t) = (vac\0(x, 0)O{y, t)|vac) c , (4.7) 

6 Actually, in Ref. the measurement was done in Fourier space (i.e. in relative momenta), but 
this difference is not significant here. For our purpose the coordinate space representation is more 
convenient. 
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Figure 4: The phase shifts in units of ty/2 vs. p/M. The open symbols are from ref. [[ll). The 
corresponding filled symbols are our measurements on the same lattices A and B, with the method 
discussed here. The crosses denote results using the FP action (D and E). The continuous curve 
corresponds to the ZZ S-matrix. 



where 

1 

0(x,t) = — ^2<r(z,t)o-(z + x,t). (4.8) 

2=0 

We omit here the O(n) structure and indices. The subscript c in eq. Q4.7|) means 
that in the 1 = channel the vacuum contribution is subtracted. 

Using the eigenvectors \n) of the transfer matrix as intermediate states one has 

Cxy(t) = J2 e ~ Ent Mx)My) , (4.9) 

n 

where 

i/j n {x) = (vac\O(x,0)\n) (4.10) 
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is the "wave function" of the corresponding state. The lowest energy states in 
eq. ( [4.9D are the 2-particle states and they will dominate at sufficiently large values 
of t. 

Note that the relative momentum 2p n of the 2-particle states is encoded not 
only in the energy E n but also in the wave function ip n (x). For the symmetric wave 
functions (/ = 0, 2 channels) one should have 

ip n (x) = A ■ cosp n (x - L/2), for R < x < L - R , (4.11) 

and similarly with sinp n (x — L/2) for the 1=1 channel. Here R is the 'interaction 
range': for a relative distance x > R the particles propagate (essentially) freely. 

One expects that the additional information obtained from the wave function 
will provide a more precise determination of p and hence of 6(p). 

4.1 Determination of the energy spectrum and wave functions. 

The rank N of the matrix C{t) in eq. @) is L/2, L/2 - 1 and L/2 + 1 in the 
/ = 0, 1, 2 channels, respectively. We assume that for t > to (with some to) no more 
than TV states contribute to C xy (t), in the sense that the total contribution of the 
states n > N is much smaller than the statistical error 5C xy (t). 

Liischer and Wolff [jll] suggested to determine the energies of the 2-particle 
states from the generalized eigenvalue problem 

C(t)v n = \ n (t,t )C(t )v n . (4.12) 

The eigenvalues of eq. fl4.12| ) are given exactly by 

A n (Mo) = e-^*-* * , (4.13) 

provided the sum in eq. (|4.9| ) is restricted to ./V terms, 1 < n < N. It is also easy to 
show that (with an appropriate normalization of v n ): 

^n{x) = C xy {to)v n {y) . (4.14) 
y 

Solving eq. ( ^4.12| ) involves an inversion of C(to), and the distortion of its small 
eigenvalues by the statistical noise is enhanced. This could affect strongly the val- 
ues and the errors of E n obtained. For to > 1 taking all N ^ L/2 states introduces 
significant instability in the result |] Because of this, we have introduced a mod- 
ification: before considering the generalized eigenvalue problem, we truncate the 

7 This equation was considered already before by Michael |28[ , in connection with a variational 
approach for evaluating the static potential in lattice gauge theory. 
8 In ref.|llj ~ L/4 states were used with to < 1. 
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correlation function to an M-dimensional subspace (M < N) spanned by the first 
M eigenvectors of C(to) (to those with the largest M eigenvalues and still stable 
against the statistical fluctuations). The generalized eigenvalue problem, eq. ( [4.12 ) 



is written then for the matrices C(t) in this reduced basis. The energies we used 
were obtained from eqs. ( |4.12| , fl~13| ). We also determined v n (x) from eq. ( P~l2D . One 
can use them as some given (nearly optimal) projectors which satisfy 

(v m , C(t)v n ) = <We- E »(*-*°) . (4.15) 

Note that due to statistical errors in C xy {t) the obtained v n (x)'s will differ from the 
true ones hence this equations will be only approximately valid. They provide a 
useful consistency check, and also a somewhat different method to determine E n . 
As an alternative way to get the phase shifts we used the momenta p n deter- 



mined from the wave function using eq. ( 4.14| ). For a given to an d M one can also 
check the self consistency of the obtained parameters by comparing C xy (t) built from 



E n and ip n (x) (cf. eqs. Q4.13J4.14|) ) with the MC result. 



4.2 The results 

We performed the calculations with the standard action for the same lattices (A,B) 
as in Ref. [|ll|] f\ In addition, we also repeated the calculations using the fixed point 
action (FP) for the 0(3) model [ p^j29[ |. The parameters of our measurements are 
summarized in Table || (here m is the inverse of the exponential correlation length) . 

The phase shifts obtained from the analysis for the standard and FP actions are 
shown in fig. || together with the results from ref. [[Tl]] . The results for the standard 
action are also given in tables |6| and [7| and for the FP action in tables |8| and ||[ The 
column <5e gives the phase shift calculated from the energy by eq. ( 4.12[ ), WF labels 



the results obtained from the wave function, eq. ( [4.14 ). The data shown correspond 



to to = 3, M = 10 for lattice D, and to = 1, M = 10 for lattice E. However, the 
results - especially for 5wf _ are quite stable against this choice. We took R ~ 3/m 
in eq. ( 4.11[ ). Taking into account that the correlation lengths used with the FP 



action were only £ ~ 3 and 6, the agreement with the analytic prediction is very 
good. Note, however, that similar suppression of lattice artifacts has been observed 
previously with this action for other observables 1 13], 29, 30]. 



Ref. [tLlP also measured an additional lattice C with mL 



21 



Acknowledgements 

We would like to thank Michael Karowski for useful discussions, and Paolo 
Butera, Paolo Rossi and Ettore Vicari for informative correspondence. 

This investigation was supported in part by the Hungarian National Science 
Fund (OTKA) (under T016233 and T019917), and also by the Schweizerischer Na- 
tionalfonds. The work of M.N. was supported by NSF grant 97-22097. 

References 

[1] K. Symanzik, Nucl. Phys. B226 (1983) 187 

For a review see M. Liischer, "Improved Lattice Gauge Theories", in Les 
Houches 1984, Proceedings, Critical Phenomena, Random Systems, Gauge The- 
ories, 359-374; and "Advanced Lattice QCD", talk given at Les Houches Sum- 
mer School in Theoretical Physics 1997, hep-lat/9802029 

[2] A. Patrascioiu and E. Seiler, Phys. Rev. Lett. 74 (1995) 1920 

[3] M. Liischer, Nucl. Phys. B135 (1978) 1 

[4] A. B. and Al. B. Zamolodchikov, Ann. Phys. 120 (1979) 253; Nucl. Phys. B133 
(1978) 525 

[5] F.A. Smirnov, 'Form factors in Completely Integrable Models of Quantum Field 
Theory', World Scientific, 1992 

[6] M. Karowski, in Field theoretical methods in particle physics, 1980, ed. W. 
Riihl, Pg. 307 

[7] M. Karowski and P. Weisz, Nucl. Phys. B139 (1978) 455 

[8] J. Balog, M. Niedermaier and T. Hauer, Phys. Lett. B386 (1996) 224 

[9] J. Balog and M. Niedermaier, Nucl. Phys. B500 (1997) 421; Phys. Rev. Lett. 
78 (1997) 4151 

[10] H. Babujian, A. Fring, M. Karowski and A. Zapletal, Nucl. Phys. B538 (1999) 
535 

[11] M. Liischer und U. Wolff, Nucl. Phys. B339 (1990) 222 

[12] J. Balog, M. Niedermaier, F. Niedermayer, A. Patrascioiu, E. Seiler and 
P. Weisz, in preparation 

[13] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B414 (1994) 785 



22 



[14] P. Hasenfratz, M. Maggiore and F. Niedermayer, Phys. Lett. B245 (1990) 522 
[15] A. Patrascioiu and E. Seiler, Phys. Letts. B445 (1998) 160 
[16] A. Pelissetto and E. Vicari, Nucl. Phys. B519 (1998) 626 

[17] M. Falcioni, G. Martinelli, M. L. Paciello, G. Parisi and B. Taglienti, Nucl. 
Phys. B225 (1983) 313 

[18] P. Butera and M. Comi, Phys. Rev. B54 (1996) 15828; ibid B50 (1994) 3052; 
ibid B58 (1998) 11552 

[19] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Phys. Letts. B402 (1997) 
141 

[20] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Nucl. Phys. B459 (1996) 
207; Nucl. Phys. Proc. Suppl. 47 (1996) 751 

[21] P. Rossi and E. Vicari, private communication 

[22] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Phys. Rev. B54 (1996) 
7301 

[23] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Phys. Phys. D54 (1996) 
1782; Nucl. Phys. Proc. Suppl. 47 (1996) 755 

[24] J. Balog, M. Niedermaier, F. Niedermayer, A. Patrascioiu, E. Seiler and 
P. Weisz, in preparation 

[25] P. Butera, private communication 

[26] J. Kim, Phys. Lett. B345 (1995) 469 

[27] S. Caracciolo and A. Pelissetto, Nucl. Phys. B455 (1995) 619 
[28] C. Michael, Nucl. Phys. B259 (1985) 58 

[29] M. Blatter, R. Burkhalter, P. Hasenfratz and F. Niedermayer, Phys. Rev. D53 
(1996) 923 

[30] S. P. Spiegel, Phys. Lett. B400 (1997) 352 



23 



